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We present a phenomenological analysis of the cos (j) and cos 2(f> asymmetries in unpolarized semi- 
inclusive deep inelastic scattering, based on the recent multidimensional data released by the COM¬ 
PASS and HERMES Collaborations. In the TMD framework, valid at relatively low transverse 
momenta, these asymmetries arise from intrinsic transverse momentum and transverse spin effects, 
and from their correlations. The role of the Cahn and Boer-Mulders effects in both azimuthal 
moments is explored up to order 1/Q. As the kinematics of the present experiments is dominated 
by the low-Q^ region, higher-twist contributions turn out to be important, affecting the results of 
our fits. 

PACS numbers: 13.88.-|-e, 13.60.-r, 13.85.Ni 


I. INTRODUCTION 

Since the early QCD investigations of hadronic hard processes, it has been recognized that azimuthal 
asymmetries in unpolarized reactions, such as Drell-Yan production and semi-inclusive deep inelastic 
scattering (SIDIS), represent an important window on the perturbative and non-perturbative aspects of 
strong interactions mm- Focusing on SIDIS processes, while at large momentum transfer Q and large 
Pt (the transverse momentum of the produced hadron), the azimuthal asymmetries are perturbatively 
generated by gluon radiation, at small Pp, Pp ^ Q, they can arise from the intrinsic motion of quarks 
[3HS] (for a review, see Ref. [7]). In the latter regime, the so-called transverse-momentum-dependent 
(TMD) factorization applies, and the SIDIS structure functions (in the current fragmentation region) can 
be written as transverse-momentum convolutions of TMD distribution and fragmentation functions [8] . 

Two azimuthal modulations appear in the unpolarized SIDIS cross section, of the type cos / and 
cos 2/, where / is the azimuthal angle of the produced hadron (measured from the scattering plane). 
These asymmetries, which have been experimentally investigated for the first time in the large region 
by the EMC and ZEUS experiments [aHii], have recently attracted a large experimental and theoretical 
attention as a potential source of information, in the small Pp region, on the so-called Boer-Mulders 
distribution function, hj-, which measures the transverse polarization asymmetry of quarks inside an 
unpolarized nucleon [I2j . 

A few years ago the azimuthal asymmetries in unpolarized SIDIS have been measured by the COM¬ 
PASS and HERMES Collaborations for positive and negative hadrons, and presented as one-dimensional 
projections, with all variables {xb, Zh,Q^, Pt) but one integrated over [T^HTS] . 

The one-dimensional data on the cos 2/ asymmetry were analyzed in Ref. m. where it was shown that 
the larger asymmetry for TT~(h~) production, compared to 7r+(/i“''), was an indication of the existence 
of a non-zero Boer-Mulders effect, in agreement with the earlier predictions of Ref. m- It was also 
pointed out that measurements at different values of were essential, in order to disentangle higher- 
twist contributions from the twist-two Boer-Mulder term. 

The HERMES and COMPASS Collaborations have recently provided multidimensional data in bins of 
a^B, z/i, and Pp for the multiplicities [mils] and for the azimuthal asymmetries EHIISI]. In principle, 
multidimensional data should offer detailed information, essential to unravel the kinematical behaviour of 
these asymmetries: for instance, the dependence is crucial to disentangle higher-twist effects. In this 
paper, we present a study of the SIDIS azimuthal moments (cos </)) and (cos 2/) in order to understand 
the role of the Cahn effect and to extract the Boer-Mulders function. We will see that, due to the present 
kinematics which is still dominated by the low-Q^ region, the higher-twist contributions are important 
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Hadron production plane 



FIG. 1: Kinematical configuration and conventions for SIDIS processes. The initial and final lepton momenta 
define the x-z plane. 


and strongly affect the results of our fits. 


II. FORMALISM 


The process we are interested in is unpolarized SIDIS: 


m + N{p) ^ i'{e) + h{Ph) + x{Px) 

The cross section of this process is expressed in terms of the invariants 

P-q P-Ph 


( 1 ) 

( 2 ) 


where q = i — t! and = —q^. 

The reference frame we adopt is the 7*-7V center-of-mass frame, with the virtual photon moving in the 
positive 25 direction (Fig. [^. We denote by Pt the transverse momentum of the produced hadron. The 
azimuthal angle of this hadron referred to the lepton scattering plane will be called (^. 

The unpolarized SIDIS differential cross section is 

dXBd„dtdPjd0 = ^{(l + (l-riVdd+2(2-rivT^FJS‘<=»0 

+ 2{l-y)F^^f^*\os{2ct>)], (3) 

where the structure functions Fjju^ P^u‘^j P^u'^'^ depend on xb, Zh,Q^, Pt- Fuu is the structure function 
which survives upon integration over tp, while F^^'^ and are associated to the cos^ and cos20 

modulations, respectively. 

If k is the momentum of the quark inside the proton, and k± its transverse component with respect 
to the "f*N axis, in the kinematical region where Pt ^ fcj, <C Q, the transverse-momentum-dependent 
(TMD) factorization is known to hold at leading twist. The structure functions can be expressed in terms 
of TMD distribution and fragmentation functions, which depend on the light-cone momentum fractions 


X = 


fc+ 
¥+ ’ 


25 = 


(4) 
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where k is the momentum of the fragmenting quark. Although the TMD factorization has not been 
rigorously proven at order 1/Q (that is, at twist three), in the following we will assume it to hold, as 
usually done in most phenomenological analyses. 

Introducing the transverse momentum pj^ of the final hadron with respect to the direction of the 
fragmenting quark, up to order k±jQ one has pj^ = —zk± and the momentum conservation reads 

Pt = zk^+p^. (5) 

At the same order we can identify the light-cone momentum fractions with the invariants Xb and 

X = XB, Z = Zh- (6) 

In the TMD factorization scheme the structure function Fuu is given by 


Fuu = ^ 

9 


d^fc_L fq{x,k±)Dq{z,p±), 


(7) 


where fq{x, k±) and Dq{z,p±) are the unpolarized TMD distribution and fragmentation function, respec¬ 
tively, for the flavor q (the sum is intended to be both over quarks and antiquarks). Note that in Eq. 0 
the transverse momentum conservation has been applied, so that p± = \Pt — zk±\. The dependence 
of all functions is omitted for simplicity. 

The structure function associated with the cost)) modulation turns out to be an order 1/Q, i.e. a 
twist-3, quantity. Neglecting the dynamical twist-3 contributions (the so-called “tilde” functions, which 
arise from quark-gluon correlations), F^^^ can be written as the sum of two terms 


pCOS 4> _ fp' 

^UU — 


iCOS (p 

UU 


F, 


Cahn 


cos (p 
UU 


BM 


( 8 ) 


with {h = Pt/\Pj 


f: 


cos p 


UU 


Cahn 


= -2E 


e^x 


2, ik±-h) 


d^k 


Q 


fq{x,k±)Dq{z,p±), 


E 


UU 


BM 


= / d^fcj. ^E /gt/p(a:,fc_L)AD^/^T(2:,pj.). 


fc_L 


(9) 

( 10 ) 


Eq. Q is the Cahn term, which accounts for the non-collinear kinematics of quarks in the elementary 
subprocess £q —)■ £'q'. Eq. (10) is the Boer-Mulders contribution, arising from the correlation between 
the transverse spin and the transverse momentum of quarks inside the unpolarized proton. In this term 
the Boer-Mulders distribution function Af^t/p couples to the Collins fragmentation function ADf^/gt. 
The relations between these functions, as defined in the present paper, and the corresponding quantities 
in the Amsterdam notation is 


k± 

^fq'^/pix, k±) ’ 

^Ph/qt(z:,P±) = ^^H^(z,p±), 


( 11 ) 

( 12 ) 


where Mp and Mh are the masses of the proton and of the final hadron, respectively. The Boer-Mulders 
effect is also responsible of the cos 2(j) modulation of the cross section, giving a leading-twist contribution 
(that is, unsuppressed in Q). It has the form 


-nCOS 2 p 

^UU 


BM 


= -E®9^ / 


PT{kB-h) + Zh [kl-2{kB-hp] 
2k±p± 


'x^fqr/p{x,k±)ADh/qt{z,P±)- 


( 13 ) 
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Concerning higher-twists, the cos 2(j) structure function has no l/Q component, but receives various l/Q^ 
contributions. Only one of these, of the Cahn type, is known and reads 


T^cos ‘2(1) 

^UU 


Cahn 


yd^fe_L 


2(fcj_ • hf 

Q^~ 


^2 

—/,(x,fc_L)n,(z,p_L). 


(14) 


Notice however that at order \/Q^ many simplifying kinematical relations do not hold any longer (for 
instance, x and z do not coincide with xb and zt), and there appear target mass corrections. Thus, 
Eq. (14) must be intended only as an approximation to the full twist-4 contribution to 

For the TMD functions we will use a factorized form, with the transverse-momentum dependence mod¬ 
eled by Gaussians. This Ansatz is supported by phenomenological analyses (see for instance Ref. [55]) and 
lattice simulations [23] . Thus, the unpolarized distribution and fragmentation functions are parametrized 
as 


e-kl/{kl) 

fq/p{Xjk±) = fq/p{x) ^^^2 ^ > ( 15 ) 

e-pl/(pi> 

Dh./q{z,p±) Df^lqi^z) / 2 \ ■ (16) 


The integrated functions, fq/p{x) and Dh/q{z), will be taken from the available fits of the world data 
(in particular, we will use the CTEQ6L set for the distribution functions [24] and the DSS set for the 
fragmentation functions [25] b The widths of the Gaussians might depend on x or z, and be different for 
different distributions: we will discuss these possibilities later on. 

For the Boer-Mulders function we use the following parametrization 


with 


Mqr/pix, k^) = A/,t/p(x) ^ ^{kl) 


^fqr/p{x) = Nq " . - x“‘>(l - xfo fq/p{x), 


(17) 


(18) 


where aq, j3q and are free parameters to be determined by the fit. The distribution so constructed 
is such that the positivity bound | A/,jt/p(a;, A:_l)| < 2 fq/p{x, k±) is automatically satisfied. Multiplying 
the two Gaussians in Eq. 0) , the Boer-Mulders function can be rewritten as 


A/,T/p(a;, k±) = A/,T/p(a;) V2e 


k± 




Tr{kl) 


(19) 


with 


(^i)BM 


(fcj) 


( 20 ) 


For the Gollins function we have a similar parametrization, namely 


AD. 


( 21 ) 


with 


AD^^qriz) = 


(7 +J)(^+^) 


zT'(l - zf Dh/q{z), 


( 22 ) 
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where 7 , <5 and Mq are free parameters. Combining the two Gaussians in (21), we get 




P± 

Me 


Q-p\l{v\)c 

7r(pi) 


(23) 


having defined 


(pi)c 


(pi)+M 2 ' 


(24) 


The Gaussian parametrization has the advantage that the integrals over the transverse momenta can 
be performed analytically. Thus, inserting the above distribution and fragmentation functions into the 
expressions for the SIDIS structure functions, we get 


2 ^~PtI{Pt) 

^UU ~ / . fq/p{Xg) Dh/q{zh) ^Ip2\ > 

q ^ ' 

i;,cos(^i _ r^Pr „2 _ f \ \ ^h{k\) e~MI{PT) 

kfju Icahn — 2 ^ 2 _^ ^B Jq/p[Xs) Ph/q\Zh) ’’"(P^) 

.cos<^ _^„Pt ST ^fqt/p{x^) ADh/qt{Zh) e-P^PPT)BM 


f: 


uu 


BM 


Pt ^ 2 

g ^ eqX_ 

. AkDLipDl 

{kl){pl) 




Me ^{P^)L 


^h{k±)BM(^PT ~ (Pt) + {p±) c{Pt)bM 


(25) 

(26) 

(27) 

(28) 


where 


and 


P^ V 

Cahn ~ ^ ^ ^ ^a^B fq/p{Xs) Dh/q{Zh) 
^ <? 


4{kir 

(29) 

7r(p2) 

(P|)2 ’ 

2 Afqr/p{Xg)ADii/gr{zh') 

BM ^Pp^^e^xs 



e Pt/{Pt)b» Zh{k\)l^{p\)l 

^{P^)L (kl){pl) ’ 



(30) 

{Pt) = {pI) + 4 (kl), 



(31) 

{Pt)bM = {p±)c + 4 {kl)BM 



(32) 


The quantities actually measured in unpolarized SIDIS experiments are the multiplicities and the az¬ 
imuthal asymmetries. The differential hadron multiplicity is defined as 

d^n^ / dVPiS \~^ dV 

dzh dP| V dy J da;^ dy dzh dP| ’ 

The deep inelastic scattering (DIS) cross section has the usual leading-order collinear expression 


dXg dy 


2 TT 


[l + (l-y)^] 

y2 


'y ^ ^qXB fqjp{Xg ) . 
(J 


(34) 
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Inserting the SIDIS cross 
the multiplicities 


section ^ integrated over (j), and Eqs. (25) and (34) into Eq. (33), we find for 


d^n'* 

dzh dP^ 

The cos (j) and cos 2(j) asymmetries 


_ Y.q elxBfq/p{x^)Dh/q{zh) e-Pil^Pj) 
'ZqeqXBfqjpixJ TT (P^) 

are defined as 

J d(j)da 


that is, in terms of the structure functions, 


co,^_ 2(2-;/)7r^f-^ 

[l + (l_y)2] 

Acos2^ ^ 

[l+(l_y)2] Fuu ■ 


(35) 

(36) 

(37) 

(38) 

(39) 


In the following we will fit simultaneously the multiplicities and the cos (j) and cos 2(f) asymmetries. 
As seen from Eq. (|35|), the data on multiplicities, which statistically dominate our dataset, constrain 


(P|.) only, which in the Gaussian model (and neglecting 1/Q^ corrections) is given by the combination 
(P|.) = + {p\), see Eq. (31). The azimuthal asymmetries, on the other hand, depend on {k\) 

and {p\) separately, offering the chance to gain information on the individual values of these two average 
intrinsic momenta. 


III. ANALYSIS OF (cos b) AND (cos 20) 

We will use the HERMES and COMPASS multidimensional data, which are provided in bins of 
XB,Zh,Q'^ and Pt for the multiplicities naiis] and for the azimuthal asymmetries [Ml mi- The HER¬ 
MES multiplicity measurements [T5] refer to identified hadron productions (7r+, 7r“, 77+, K~) off pro¬ 
ton and deuteron targets and cover the kinematical region of values between 1 and 10 GeV^ and 
0.023 < Xg < 0.6, while the COMPASS multiplicity data [19] refer to unidentified charged hadron pro¬ 
duction (h+ and h~) off a deuteron target (®LiD), cover the region 0.0045 < < 0.12 and are binned in 

a region similar to that of the HERMES experiment. In Ref. [Mjj the HERMES collaboration released 
azimuthal asymmetries for unpolarized SIDIS identified hadron production off proton and deuteron tar¬ 
gets. Multidimensional data on the and asymmetries are provided, binning the data in the 

kinematical variables Xb, y, zu and Pt- The kinematical cuts used for their analysis are the following 

0.023 <xb< 0.6, 0.2 < Zh < 1.0, 0.05 GeV^ < Pt < 1.3 GeV^ , (40) 

with the additional constraints 

> 1.0 GeV^, > 10.0 GeV^ , 0.2 < y < 0.85 (41) 

binned according to Table H of Ref. m- One dimensional projections of the azimuthal moments are also 
provided, integrated on the kinematical regions presented in Table HI of Ref. [Mj • 

The COMPASS collaboration azimuthal asymmetries for unpolarized SIDIS charged hadron-production 
off a deuteron target are presented in Ref. m- In their analysis, both one-dimensional and multidimen¬ 
sional versions of and were extracted, binning the data in the kinematical variables xb, Zh 

and Pt- The kinematical cuts used for their analysis are the following 

0.003 <XB < 0.13, 0.2 <zh< 0.85 , 0.1 GeV^ < Pt < 1.0 GeV^ , 


(42) 
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with the additional constraints 


g2 > 1.0 GeV^ , > 25.0 GeV^ , 0.2 < y < 0.9. 


(43) 


In the following, we are going to study the multidimensional SIDIS azimuthal moments (cos (/)} and 
(cos 2(/>} with the aim of understanding the role of the Cahn and Boer-Mulders effects. We expect multi¬ 
dimensional data to be useful to explore the kinematical behavior of these asymmetries, and their explicit 
dependence to be of help in exploring the interplay between leading-twist and higher-twist contribu¬ 
tions. We perform our analysis by including contributions up to order 0(1/Q); the role of dynamical 
twist-three contributions will be briefly discussed at the end of the Section. 

Up to order 0(1/Q), (cost/)) receives contributions from the Cahn and the Boer-Mulders effect, which 
appear at subleading twist in Eqs. ([Q) and (10), while (cos2^) is proportional to the sole Boer- 

Mulders effect, which instead appears at leMing twist, Eq. (13). Both asymmetries contain, at denomi¬ 
nator, the contribution of Fuu, defined in Eq. ([^. 

Fjjij and the Cahn contribution to (cos (f>) involves only the unpolarized TMD distribution and fragmen¬ 
tation functions fq/p{x,k±) and Di^^q{z,p±). These functions have been recently extracted in Ref. [55], 
from a best fit of HERMES and COMPASS multidimensional multiplicity data. In this analysis, a 
Gaussian model was used for the kj_ and p± dependence as in Eqs. (15) and (16). 


If we use the values of {k^) and {p\) obtained there we find a very large Cahn contribution, of the order 
of 50% which largely overshoots the data. This is not surprising. Eq. (26) shows in fact that E'))y'^|cahn 
is proportional to {k'j_), which was found to be rather large, about 0.6 GeV^. It is quite unlikely that any 
reasonable Boer-Mulders contribution could cancel this huge Cahn term so as to reproduce the observed 
(cos (/)) asymmetry, which does not exceed 10%. 

However, in Ref. |26j it was o bserved that, since the multiplicities are sensitive only to the combination 
(P|.) = z‘l{k‘j_) + {pY), Eq. (25), they cannot distinguish (fc^) from (p^). Instead, the azimuthal asymme¬ 
tries involve (A:^) and (p^.) separately, and are sensitive to a z^-dependent (p5_), see Eqs. (26)-(30). For 
example, if one takes {k'j_) as a free constant parameter, (fc^) = C, and allows (p^) to have a quadratic 
Zh dependence of the form 


(pi) — ^ + Bzl , 

where A and B are two additional constant parameters, one finds 

(p|.) = A + {B + C)z1 , 


(44) 


(45) 


which is of the same form of Eq. (31). A fit of the multiplicities using this functional form would lead to 
the same results, but would allow for a different interpretation of the extracted parameters in terms of 
{k\) and {p\). Instead if, in addition, we fit also (coscj)) and (cos2^) data, we acquire some degree of 
sensitivity to the individual values of A, B and C. Here, in principle, C can be small and lead to a much 
smaller Cahn contribution. 

In this paper we explore this configuration. We perform a global best fit which includes the multi¬ 
plicities, the cos(/) asymmetry and the cos2(/) asymmetry. Working up to order 1/Q, these asymmetries 
read: 


^COS 4> _ ^cos c, 


\ Cahn 


^COS cj) 


BM 


acos2(;)> _ acos2(;)>| 

~ Ibm 


As both COMPASS and HERMES data on (cos (/>) and (cos 2(j)) are restricted to a narrow x range, they 
do not allow us to determine the precise ^-dependence of the Boer-Mulders function. Thus we take 
A/, t/p to be simply proportional to fq/p, by setting Uq = (3q = 0 va. Eq. (11^. For the Collins function, 
we distinguish a favored and a disfavored component, and we fix their parameters to the values obtained 
in a recent fit of the Collins asymmetries in SIDIS and e“'"e“ annihilation [27] : 


Ng, = 0.49, iv£f = -I.OO, 
7 = 1.06, 5 = 0.07, 

= 1.50 GeV^. 


(46) 














FIG. 2: Best fit curves for (costj)} obtained by fitting COMPASS data on multiplicities, {coscp} and (cos 20) . The 
Cahn effect in (cos 20) has been set to zero. 


To parameterize the Boer-Mulders and Collins functions we need to input the unpolarized fq/p{x) and 
Dq/p{z), see Eqs. (18) and (23). Consistently with our previous choice, Eqs. (IT^ and (T^, we will use 
the collinear CTEQBL distribution functions [23] and DSS fragmentation functions [25 . As mentioned 
before, (fc^) = C and {p\) = A + Bz^, with A, B and C free parameters to be determined by the fit. 

It is known that the COMPASS multiplicities should be corrected by a large normalization factor: in 
fact, issues in that analysis were detected, which can affect the overall x, y, z normalization of multiplicities 
up to 40%, as pointed out in Ref. |2H]- Lacking for the moment the corrected data, in the present paper we 
apply the same multiplicative normalization factor as obtained in [26] . Ny = 1.06 — 0.43?/. This correction 
was found to improve considerably the fit of the COMPASS multiplicities. 

The kinematical range explored by the two experiments is further restricted in order to make sure 
that our description, based on TMD factorization, can safely be applied. To avoid contaminations 
from exclusive hadronic production processes and large 0 resummation effects [29] we select data with 
z/j < 0.65 for COMPASS and with Zh < 0.69 for HERMES. The lowest cut in is chosen accordingly 
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FIG. 3: Best fit curves for (cos 20) obtained by fitting COMPASS data on multiplicities, (cos0) and (cos 20) . The 
Cahn effect in (cos 20) has been set to zero. 


to the minimum in the CTEQ6L analysis, > 1.69 GeV^, which amounts to excluding the lowest x 
bins. Finally, we select 0.2 < Pt < 0.9 GeV, following Ref. pSj . 

The results of the fit for the azimuthal moments are shown in Figs. [^and[^ for the COMPASS data, 
and in Figs. I^and for the HERMES data. The description of the multiplicities is practically unchanged 
compared to Ref. [26j . therefore we do not show the plots here. The values of of the parameters 

are listed in Tables |T] and IE 

The asymmetry data (especially (cos 0)) drive the transverse momentum of quarks to a very small 
value, (A:^) ~ 0.03 — 0.04 GeV^, which means that the transverse momentum of the produced hadron is 
largely due to transverse motion effects in the fragmentation process. The difference between positive and 
negative hadrons is found to be small (or even negligible), and the agreement with the data on (cos0) 
worsens as Zh grows and decreases. 

One may wonder whether a more flexible model, including flavor dependent Gaussian widths, could 
modify our results. It would indeed be interesting to determine whether the available SIDIS data signal 
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TABLE I: Minimal and parameters, for a fit on COMPASS multiplicities, (cos (j)) and (cos 2(/!)) in which the 
Cahn effect in (cos 2(j)) has been set to zero. Parameter errors correspond to a 2a confidence level. 


Cuts 



Parameters 


Zh < 0.65 

iXpt'jmult ~ 

3.43 

A = 0.200 ±0.002 

MIm = 0.09 ± 0.45 

> 1.69 

{Xpt')cos(4>) 

1.17 

B = 0.571 ±0.018 

Nd = -1.00 ± 1.95 

0.2 < Pt < 0.9 

(Xpt)cos(20) “ 

= 1.02 

C = 0.031 ±0.006 

Nu = -0.45 ± 0.26 

Minimal parameters, for a fit on HERMES multiplicities, (cos 0) and (cos 20) i 

in (cos 20) has been set to zero. Parameter errors correspond to a 2a confidence level. 

Cuts 

x^ 


Parameters 


Zh < 0.69 

{,\pt')mult — 

1.70 

A = 0.126 ± 0.004 

MIm = 0.10 ±0.20 

> 1.69 

{Xpt)cos(4>) — 

: 2.39 

5 = 0.506 ±0.045 

Nd = -1.00 ±0.20 

0.2 < Pt < 0.9 

{.Xpt'}cos(2(p) ~ 

= 2.13 

C = 0.037 ±0.004 

Nu = -0.49 ±0.15 


any flavor dependence in the unpolarized TMDs. This was already considered in the analysis of multi¬ 
dimensional multiplicities of Ref. |26j . where it was noted that flavor dependence did not improve the 
description of the data significantly. Here we have performed several fits of the asymmetries allowing for 
flavor dependent parameters, but we have found that the overall picture does not improve. Furthermore 
including flavor dependence generates an over-parameterization, given the precision of present data, and 
consequently results in largely unconstrained fit solutions. 

From Tables and one sees that the presence of the Boer-Mulders function is rather marginal: 
Afdt/p is very uncertain and compatible with zero, whereas A/„t /p is slightly more constrained, but very 
small. 

The reason of this result is that our selection of multidimensional data on cos 20, which cuts out a large 
portion of data corresponding to small values, turns out to be compatible with a zero asymmetry. 
Notice however that in the small and large z region the asymmetry is instead quite sizable. This can 
be seen also by considering the previous, one-dimensional data, where all variables but one are integrated 
over [I3HI1]- This means that the integrated asymmetries are mainly driven by small and large z 
events, which could be affected by relevant higher-twist contributions. The importance of the 
Cahn term, Eq. (291, was indeed pointed out in Ref. [16], but one should not forget that this term is 


only a part of the overall twist-four contribution, which is not explicitly known (at this order there are 
also target-mass effects, and the identification oi xb and zu with the light-cone ratios x and z is no more 
valid). 

By inspecting the recent release of azimuthal asymmetries by HERMES [20] and COMPASS [21], the 
difference between the sizes of the negative and positive hadron production asymmetries is evident in the 
one-dimensional data (see Figs. 14-16 of Ref. [20] and in Figs. 14, 15 of Ref. m) ; moreover, in this case 
both and are consistently different from zero. However, these same features are not readily 

visible in the multidimensional data sets. 

At this stage we would like to understand whether we can assign a physical meaning to the extracted 
parameters. As a matter of fact, the values of the average momentum (fc^) extracted from our fit 
is quite different (one order of magnitude smaller) from that previously extracted in several different 
fits |26L 1301131j . This small value is mainly driven by the (cos 0) asymmetry, which is entirely twist-three. 
It is clear that one should make sure to evaluate this asymmetry correctly in order to interpret the 
extracted parameter as the average momentum (/c^). 

The Cahn and Boer-Mulders terms in (cos0) can, in fact, be supplemented by other dynamical 1/Q 
contributions [32] . Taking quark-gluon correlations into account Eq. acquires, for instance, an addi- 
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<y> 


A 


<2cos((|))> TC"^ from P 



Pj, (GeV) 


► 


0.06 


0.10 


0.19 


0.36 


FIG. 4: Best fit curves for (cos0) obtained by fitting HERMES data on multiplicities, {coscj)) and (cos 20). The 
Cahn effect in (cos 20) has been set to zero. Here we show only one bin in 2 , as an example, with {z) = 0.45 


tional term containing a new distribution function, /g, and becomes 


F, 


uu 


Cahn 


= -2^elx / 


Q 


[fqix, k±) + fqix, kj_)]Dq{z,p±). 


(47) 


Notice that this new term cannot be separated from the other. A similar correction, with another 


distribution function hq, applies to the Boer-Mulders component (101, and the fragmentation functions 


should also be modified for quark-gluon correlations. 

Whereas it is generally believed that the “tilde” 0{1/Q) contributions are negligible, it is possible that 
in the kinematical regions presently explored this might not be the case. Possible cancellations among 
these terms could affect the extraction of the intrinsic momenta parameters. 

In order to estimate the impact of the dynamical higher-twist contributions to the Cahn effect, we can 
simply assume tilde functions to be proportional to the non-tilde ones. Given the restricted kinematical 
ranges of the data, this is not a very limiting assumption. Thus the presence of fq is effectively simulated 
by an extra normalization constant Ncahn in front of the (cos (jj) Cahn term, that is: 


^COS0 ^ kl™^1cahn + 


IBM 
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FIG. 5: Best fit curves for (cos2(^) obtained by fitting HERMES data on multiplicities, {coscp} and (cos 20) . The 
Cahn effect in (cos 20) has been set to zero. Here we show only one bin in z, as an example, with (z) = 0.45 


Since (cos0) is dominant in our fit, one sees from Eq. (26) that the effect of Afcahn is compensated by a 
readjustment of (fc^) (which determines not only the width of the Gaussian distributions, but also the size 
of the asymmetries). Therefore, identical fits are obtained by allowing fVcahn to be smaller than unity, 
and proportionally increasing For instance, setting fVcahn = 0.5, that is assuming that dynamical 

twist-3 terms reduce by 50% the Cahn term, one gets {k‘j_) = 0.06 GeV^ (twice the value obtained in our 
main fit). An even larger cancellation, that is a smaller A"cahn ~ 0.1 coefficient, would deliver a value of 
(/c^) similar to those extracted in previous analyses P?1 IHT] . 

We conclude that in the present kinematics the structure and the magnitude of the higher-twist terms, 
which are not fully under control, are crucial for determining (fc^). 
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IV. CONCLUSIONS AND PERSPECTIVES 

In the TMD framework, the (cos </>) and (cos 2^) asymmetries are sensitive to the transverse momentum 
of quarks inside the target and in the fragmentation process. In the Gaussian model of quark distributions, 
the widths (fc^) and {p\) also determine the size of the asymmetries. Adopting an 0{1/Q) scheme, 
which attributes (cos (/)) both to Cahn and Boer-Mulders effects at order 1/Q, and (cos2(^)to the Boer- 
Mulders effect at leading twist, and ignoring twist-3 dynamical contributions (arising from quark-gluon 
correlations), our analysis shows that the recent COMPASS and HERMES multidimensional data can 
be reproduced by a very small value of {k\), namely 0.03-0.04 GeV^. Within this picture, this means 
that most of the transverse momentum of the outgoing hadron is due to the fragmentation, which must 
be described by a function with a z-dependent width. This result, mainly driven by (cos^), could be 
modified by the presence of further twist-3 terms, which might not be negligible due to the relevance of 
the small-Q^ region in the present measurements. 

A somehow disappointing output of our fits is the indeterminacy on the extraction of the Boer-Mulders 
function, which seems to play a minor role in the asymmetries. This is seen in particular from (cos 2(^), 
which is entirely determined by the Boer-Mulders contribution but appears to be, within large errors, 
compatible with zero. 

On the other hand, the integrated (cos 2^) data [ID] show a non vanishing asymmetry, especially when 
plotted against z. The asymmetry is slightly negative for 7r+ and positive for 7r“, as expected from the 
Boer-Mulders effect m- Also the integrated data on (cos (p) show a different asymmetry for 7r+ and 
TT~: this indicates a flavor dependence which can only be achieved with a non-zero Boer-Mulders effect 
since, within a flavor-independent Gaussian model with factorized x and k± dependences, the Cahn 
effect is flavor blind and can only generate identical contributions for positively or negatively charged 
pions. However, the signs of the u and d Boer-Mulders functions required for a successful description of 
(cos 2(j)) appear to be incompatible with those required to generate the appropriate difference between 7r+ 
and TT~ in the {cosp) azimuthal moment. As we mentioned, a more refined model with flavor dependent 
Gaussian widths is not helpful, given the precision of the current experimental data. 

One should not forget about the existence of other higher-twist effects that could combine with the 
Boer-Mulders term and alter the simple picture considered here. In order to disentangle these contribu¬ 
tions, it might be useful to integrate the asymmetry data on restricted kinematical ranges, so as to avoid 
the low-Q^ region and meet the requirements of TMD factorization. Analyzing properly integrated data 
could help to clarify the origin of azimuthal asymmetries and possibly to get more information on the 
Boer-Mulders function. Work along these lines is in progress. 

It would also be interesting to investigate how SIDIS azimuthal modulations can be affected by gluon 
radiations. Following Ref. |6] one can actually compute the perturbative corrections originating from 
gluon radiation at order for the numerator and denominator of the azimuthal asymmetries. Indeed, 
in the limit of small qx (where our analysis applies) they are affected by strong divergences, generated 
by soft and collinear gluon radiation. One might expect that these divergences cancel out in the ratio 
when building the azimuthal asymmetry, but this only happens when the divergences appearing in the 
numerator and in the denominator are of the same origin and have the same structure. In fact, in a more 
recent paper |DD], it is explicitly shown that, for Drell-Yan scattering processes, the (cos 2^) azimuthal 
modulation of the cross section shares with the azimuthal independent (integrated) term the same qx 
logarithmic behaviour of the asymptotic cross section, proportional to /q^) \n{Q^/q'j.). In this case 

the same resummation techniques which are known to work for the integrated cross section (which 
appear at denominator) can be applied to the azimuthal modulation appearing at numerator. However, 
they point out that this does not happen for the (cos (p) azimuthal modulation, which does not exhibit 
the usual diverging logarithmic term, but is simply proportional to (Q^/qx)- In this case, the usual 
resummation scheme techniques cannot be applied and new strategies needs to be devised. To the 
best of our knowledge, this has not been explicitly studied for SIDIS processes, where applying usual 
resummation schemes and conventional matching recipes is much more problematic, even in the simplest 
case of integrated, unpolarized cross sections [51] . 
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